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EDM is a well-established hole maehining option with various advantages 
due to non-eontaet eharaeteristies of the proeess. However, knowledge about 
the proeess is not enough for its more improvements. Exprimenal studies are 
eostly and time eonsuming beeause of the eomplex nature of proeess. 
Therefore, proeess modeling is a good alternative to reduee the experimental 
expense related to the teehnology. This paper studys EDM proeess through 
mathematieal model, whieh ineludes the preeise insight into the interaetive 
behavior of EDM system. The ignition, diseharge and reeovery phases of the 
model have been developed through MATLABs time domain analysis. 
Simulation result shows good agreement with expeeted profile of EDM 
spark. To verify the model, simulated material removal rates (MRRs) from 
series of simulation are eompared with the experimental ones reported by 
previous researeher. Ability of the model to prediet the dynamie behavior 
profile of the EDM system is sueeessfully eonfirmed by low average 
pereentage error in predieting MRR. 
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I. INTRODUCTION 

EDM is a non-traditional machining process to remove material from the workpiece using a thread 
of electrical discharges between the electrode and the workpiece called gap [1]. non-contact characteristics of 
EDM makes it a valuable technique for variety of hole manufacturing approaches [2]. However, scientific 
knowledge about the process is still insufficient and it is the main obstacle for its more improvements. 
Experimental trials are challenging due to the highly stochastic and complex nature of the process caused by 
some process complications such as adhesion, short-circuiting and cavitations that increase the machining 
time and make the process become unstable [3]. 

Many EDM researchers have attempted innovative ideas such as ultrasonic vibrations and flushing 
effect with the purpose to solve the process complications and improve the system stability. Mahardika [4] 
found that using ultrasonic vibrations remove adhesion and short circuit during machining procedure. 
Shabgard [5] revealed that EDM of Ti-6A1-4V by using ultrasonic vibrations of copper electrode enhace 
MRR. He also claimed that crack density and tool wear ratio (TWR) reduce at finishing regimes while recast 
layer, cracks density, and TWR increase at roughing regimes. Although, ultrasonic vibrations help to promote 
productivity and feasibility in certain conditions, however there is still a long way to go for use in industry. 

Flushing transfers dielectric fluid into the pipe electrode, so debris carries away and the insulating 
characteristics of the dielectric are preserved, it is an easy way to improve efficiency of EDM hole drilling 
[6]. Pressure and suction flushing techniques have the capability of eliminating the spark eroded particles, but 
they can only be used in particular applications. Ebisu [7] investigated the effect of jet flushing on the debris 
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stagnation during the cutting corner shape in Ist-cut wire EDM. He showed that applying jet flushing 
decreases debris stagnation in the gap only temporarily after the corner. 

In order to overcome the difficulty aspects related to real machining environment and its effect on 
system stability, modeling and simulation is an alternative way to understand the mechanism underlying the 
EDM procedure. Although gap spark profiles play an important role in the analysis of the EDM process, but 
very few mathematical models have been developed to identify the profile of EDM spark during machining 
procedure. Minhat et al. [8] presented the mathematical model of EDM pulses based on the initial, ignition 
and discharge phases. However, the model failed to determine the EDM profiles correctly. Moreover, 
mathematical equation did not match the profiles in all phases. 

This paper proposes mathematical model to predict the dynamic behavior of the EDM spark, similar 
to the ideal one. The mathematical explanation of the model is located in section 2. Ignition, discharge and 
recovery phases have been considered properly. In section 3, MATLAB software is used to simulate 
machining procedure. Section 4 includes discussion about the validity of the simulated model using series of 
experimental data to calculate MRRs and compare with the experimental results from previous researcher. 
Finally, conclusion is given in Section 5. 


2. MODEL DESCRIPTION 

In this section, mathematical model during ignition, discharge and recovery phases is developed 
separately. The proposed model shown in Figure 1 consisted of pulse power generator and EDM spark. Pulse 
power generator, in turn, is made of an ac source, a transformer, rectifier and a capacitor as filter followed by 
a transistorized switching circuit as pulse generator. All components are supposed to be ideal in order to 
reach simple and clear insight into the model behavior. 

Based on the model, AC source voltage is applied to provide required value of DC input voltage 
Vin through the bridge diode and transformer with the turns ratio of ni:n 2 where ni and n 2 are the numbers of 
primary and secondary windings, respectively. Pulse generator comprised of a MOSFET switch Si and a 
resistor Ri. Switch Si is controlled by low power pulses of P^i and used to control current from capacitor C 
to the gap. As can be seen in Figure 1, the EDM spark which is the gap model between electrode and 
workpiece consists of Rs, Rig, Rdis, LcHs and switch S 2 driven by Pw 2 . Three defined phases of EDM pulses with 
related mathematical equations are explained in following. 


2.1. Ignition phase 

Based on the profile illustrated in Figure 2, ignition phase is occurred in the time interval from ti to 
t 2 which is called delay time td. Open gap voltage Voc provided a severe electric field between electrode and 
workpiece. An ionization path created through the dielectric and flow of the current igap is interrupted. Gap 
voltage Vgap in this phase is equal to Voc- Small delay time results larger spark time so more energy enters into 
the workpiece [9], . Equivalent model of EDM system in this phase is obtained from Figure 1 when switch Si 
is close and switch S 2 is open. Related equations of gap current igap and gap voltage Vgap obtained as 
following. As S 2 is considered open so, 

(1) 


( 2 ) 


substituting (2) in (1), 


^gap ^Rig 


using Ohm’s law. 


V 

gap 

gap r> 

ig 


R. 


V 

_ gap 


( 3 ) 


Switch Si is on, thus applying voltage driver between resistors Rig and Rs gives: 
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V =V (■ 

sap in V 






By substituting (3) in (4), 


V =V ( 

sap in V 


V„ 


^gap ^s^gap 


Simplifying both sides of (5) gives, 


y. -y 

in gap 

R 


= i 


gap 


( 4 ) 


( 5 ) 


( 6 ) 


From (6) ean be seen that small differenee between and Vgap makes the model more elose to the 
ideal value of igap in this phase whieh is zero. So by applying (7), gap voltage reaehes to its maximum value 
ealled open gap voltage Voc • 


Then, 


y «y 

gap in 


hap “ 0 


( 7 ) 

( 8 ) 


Defieieney of this assumption in [8], led to a mismateh of its mathematieal expression with the gap 
eurrent profile in this phase. 


PVilsc Power OcDcralor 

< -► 

Powisr Supply 



Figure 1. The EDM model in ignition phase 


2.2. Discharge phase 

The diseharge phase is oeeurred in the time interval from t 2 to tg whieh is ealled tdh as shown in 
Figure 2. During this phase, the isolating effeet of the dieleetrie breaks down, eurrent starts to flow while the 
voltage falls [10]. The spark is formed and maehining eontinued to reaeh a peak gap eurrent of Ig and a 
diseharge voltage of Vdis- In this phase both switehes of Si and S 2 from equivalent EDM model in Figure 1 are 
eonsidered on. Ldis and Rdis are eonneeted together in series whieh both are eonneeted to resistanee Rig in 
parallel. Gap eurrent igap and Gap voltage Vgap ean be obtained as follows. 


^gap hj,, + \ 


( 9 ) 
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y _ / D , r ^^Rdis 

* gap ^dis ^dis dt 


modifying (10) gives, 


di. 


^dis 


dt — . 

^dis ^gap ~^R(iis ^dis 


Integrating both sides of the (11), 


'3 = r 

'3 Jo 


^'Rdis dij 


Rdis 


By applying assumption. 


^ ^gap ^Rai,^dis 


di — —^ 

^dis ^dis 


So, — can be explained as follows: 




dz 

z 


?3 _ C^Rdis 

^2 “Jo 

term (15) can be explained as follows, 
_^\t3 =inA/ -i R 

Ldis \t 2 %is-^disJo 

using limits, 

_ ^dis(h~h) _ -^^ ^gap~^Rdis^dis ^ 

^dis ^gap 

taking antilog of both sides of equation (17), 


Rdis(^2 ^2^ 

y e =V -i R 

^ gap^ ^ gap ''Rais dis 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


the current can be explained as follows: 

“ R,i, ^ ’ 

also. 


(19) 


( 20 ) 


By inserting equations (19) and (20) in equation (9), the current gap in this phase can be found as follows: 




( 21 ) 


applying the Kirchhoff law. 




( 22 ) 
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substituting (21) in (22), 

^dis (^3 “^2 ^ 

^ gaplR,^,y^ ^ 

Vgap is obtained as follow: 

^ ^dis _w 

^ gap ~ Rdis ^ in 

R,,R,{\-e ^ ) + RsRdis+Vdis 


(23) 


(24) 


From (21) and (24), it is clear that, final equation of gap current and gap voltage in this phase is 
quite different from that of presented in [8]. Expression of gap current and gap voltage in [8] can not describe 
the dynamic behavior of system accurately. 


2.3. Recovery phase 

The recovery phase is happened during the time interval from ts to t 4 which is called tree as shown in 
Figure 2. The flow of current is stopped and desired insulating electric properties of the dielectric fluid are 
recovered [11]. The schematic circuit of the EDM model in this phase is obtained from Figure 1 when switch 
Si is off and no current goes through Ri and Rs. So Vgap and igap are equal to zero. This phase is totally missed 
in the model presented in [8]. 

Since the switch is open, so: 


y 

gap 

i 

gap 


= 0 
= 0 


(25) 



Figure 2. Profile of Vgap and igap during one spark cycle over time 


3. SIMULATION 

In this section, the diagram of the EDM system with ac to dc power supply and transistorized 
switching circuit as pulse generator designed in MATEAB and schematic diagram is shown in Figure 3. 
According to the estimation obtained in (7), transformer reduces the source voltage of 250 y to the input 
voltage of Vin as near as open gap voltage Voc of 160V. Figure 4 shows simulation results of gap voltage and 
gap current for a selected machining process. Simulated data are chosen from previous experimental tests 
presented by A. Yahya. [12]. In order to optimize discharge conditions and based on the experimental test 


Int J Pow Elec & Dri Syst, Vol. 10, No. 2, June 2019 : 874 - 881 



























Int J Pow Elec & Dri Syst 


ISSN: 2088-8694 


□ 879 


[12], delay time td is set to Ifis which is insignificant compared with the discharge time tdh- It is clearly seen 
in Figure 4 where simulation results are quite correlated with desirable profile of Vgap and igap in Figure 2. 




(1 

0.1 


0.12^3 0.12^^55 0.12% OLllVeS 0.12V7 0.12^75 0.]2Wi 0.]2^HfJ 0.l2W 

] illK (si 


Figure 4. Spark profiles for Ig =12.5 A and Fs =17.24 KHZ 


4. RESULT AND DISCUSSION 

As shown in Figure 2, the profile of EDM spark during one cycle consists of the initial phase 
occurred during the time interval between ti and ^ 2 , followed by the discharge phase from t 2 to ts and the last 
phase which is recovery from ts to td. Conforming to the equation (7), small difference between input voltage 
and gap voltage is selected to get the model close to the ideal profile through ignition phase. Noteworthy 
point to mention is that to better evaluate the procedure, proposed model did not considered noise during 
EDM process which results from the stochastic nature of the EDM spark. MATFAB software is used to 
develop the complete model of EDM system. Simulation results in Figure 4, as gap voltage and gap current, 
are absolutely similar to the modeling profile from Figure 2. 

To verify the simulated model, predicted MRRs from series of simulations are compared with series 
of experimental MRRs carried out by A. Yahya [12] using steel workpiece and copper electrode. Predicted 
MRR is determined by inserting simulated data into equation (26) obtained by the same researcher in [13]. 


MRR=aVJd,,F, 



ft ^ 

3 

N 

2 

(t ^ 

‘'dis 


3.52x10“’ 

‘'dis 

-1.33x10“ 

+ 1.25x10“ 

+ 1.53 


J 




\^d J 



(26) 


where a is material properties factor. For the present study, a eonsidered to be equal to the one seleeted by 
[14] i.e., ct]=2xl(T^V/"f Equation (26) is valid for ^jj^up to 400 iis and td equal to 2 iis. All data in current 
research conform to this range of validations. 

Table 1 presents the simulated (Predicted) MRR and the experimental (Actual) MRR during several 
peak gap currents Ig , discharge times tdis, recovery time tree and spark frequencies Fs. The last column of this 
table shows predicted error which is a comparison between the experimental and the simulated MRR under 
identical conditions. The average simulated error is below of 5.14 %. It seen that the simulated model has 
ability to predict the MRR with acceptable error. 
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Table 1. Compparision between experimental (Actual) and simulated (Predicted) MRR 


Process 



Kec(j^) 

FJkHz) 

MRR {mnd / min) 

Actual Predicted 

Error (%) 

1 

4 

2 

4 

125 

4 

4.48 

12.00 

2 

4 

3 

4 

111.1 

6 

6.05 

0.83 

3 

4 

4 

4 

100 

8 

7.36 

8.00 

4 

4 

6 

4 

83.33 

10 

9.34 

6.60 

5 

4 

12 

4 

55.55 

13 

12.62 

2.92 

6 

4 

25 

4 

32.25 

15 

15.56 

3.73 

7 

4 

50 

6 

17.24 

17 

18.14 

6.71 

8 

4 

100 

12 

8.77 

19 

19.52 

2.74 

9 

4 

200 

25 

4.41 

13 

14.00 

7.69 

10 

4 

400 

50 

2.21 

12 

11.77 

1.92 

11 

6 

2 

4 

125 

7 

6.87 

1.86 

12 

6 

3 

4 

111.1 

9 

9.01 

0.11 

13 

6 

4 

4 

100 

11 

11.07 

0.64 

14 

6 

6 

4 

83.33 

12 

13.94 

16.17 

15 

6 

12 

4 

55.55 

19 

18.97 

0.16 

16 

6 

25 

4 

32.25 

23 

23.49 

2.13 

17 

6 

50 

6 

17.24 

26 

26.44 

1.69 

18 

6 

100 

12 

8.77 

21 

21.97 

4.62 

19 

6 

200 

25 

4.41 

23 

21.44 

6.78 

20 

6 

400 

50 

2.21 

19 

18.07 

4.89 

21 

8.5 

3 

4 

111.1 

11 

11.83 

7.55 

22 

8.5 

4 

4 

100 

16 

15.44 

3.50 

23 

8.5 

6 

4 

83.33 

21 

19.39 

7.67 

24 

8.5 

12 

4 

55.55 

23 

24.36 

5.91 

25 

8.5 

25 

4 

32.25 

31 

30.87 

0.42 

26 

8.5 

50 

6 

17.24 

36 

37.58 

4.39 

27 

8.5 

100 

12 

8.77 

38 

40.46 

6.47 

28 

8.5 

200 

25 

4.41 

33 

36.61 

10.94 

29 

12.5 

3 

4 

111.1 

16 

14.54 

9.13 

30 

12.5 

4 

4 

100 

20 

21.17 

5.85 


5. CONCLUSION 

In this paper, a time domain mathematical model of EDM system has been developed. The whole 
model is simulated in accordance to the EDM conditions including ignition, discharge and recovery phases. . 
MATLAB simulation result is quite correlated with desirable spark profiles. Validity of the simulated model 
is carried out by comparing MRR from the previous researcher’s experimental results. It is found that, well- 
designed model for EDM system can easily provide the possibility to predict dynamic behavior of pulse 
profiles by eliminating complications related to the stochastic nature of EDM process in experimental trial. 
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